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Abstract 

Many logic programming based approaches can be used 
to describe and solve combinatorial search problems. 
On the one hand there are definite programs and con- 
straint logic programs that compute a solution as an 
answer substitution to a query containing the variables 
of the constraint satisfaction problem. On the other 
hand there are approaches based on stable model se- 
mantics, abduction, and first-order logic model gener- 
ation that compute solutions as models of some the- 
ory. This paper compares these different approaches 
from point of view of knowledge representation (how 
declarative are the programs) and from point of view 
of performance (how good are they at solving typical 
problems). 



Introduction 

Consistency techniques are widely used for solving finite 
constraint satisfaction problems (CSP). These tech- 
niques have been integrated in logic programming, re- 
sulting in finite domain constraint logic programming 
(Van Hentenryck 1989). In this paradigm, a program 
typically creates a data structure holding the variables 
of the CSP to be solved, sets up the constraints and 
uses a labelling technique to assign values to the vari- 
ables. The constraint solver uses consistency techniques 
(Tsang 1993) to prune the search. Solutions are given 
by an answer substitution to a goal. This leads to a 
rather procedural programming style. Moreover, the 
problem description is not very declarative because the 
mapping between domain variables and their value has 
an indirect representation in a term structure. A more 
detailed discussion for representing and solving CSP in 
different logical systems, following this approach, can 
be found in (Mackworth 1992). 

An alternative representation of a CSP can be given 
in a first-order logic setting, where instead of using vari- 
ables, one can define the problem in a more natural way 
by using constant and function symbols and specify the 
constraints as logic formulae. A solution is then given 
by a model of the theory, and in particular by the inter- 
pretation of the function symbols. We argue that repre- 
senting a problem in this way tends to be more declar- 
ative. To achieve a similar declarative representation 



in a logic programming system, the functions should be 
replaced with predicates and then an answer of a CSP is 
given by a table of facts. Abduction (Kakas, Kowalski, 
& Toni 1992) is one framework which allows such rea- 
soning. The relation between the variables of the CSP 
and their values is declared as an abducible and solu- 
tion of the problem is given by a set of abduced atoms. 
More recently, a logic programming paradigm based on 
stable model semantics (Gelfond & Lifschitz 1988) has 
emerged and Niemela (Niemela 1999) proposes it as a 
constraint solving paradigm. In this approach a solu- 
tion is given by a stable model of the program. 

In this paper, we use two typical CSP problems to 
compare the merits of the various approaches. One 
is graph coloring where the domain size is constant 
but the number of constraints increases with increasing 
problem size; the other is the n-queens problem where 
both the domain size and the number of constraints in- 
creases with increasing problem size. For both problems 
we show the representation in the various approaches^, 
comment on the declarativity of the representation and 
briefly describe the basic principles behind the imple- 
mentation. We also compare the performance of differ- 
ent systems and show how it scales by increasing the 
problem size. 

First, we give some basic notions about constraint 
satisfaction and logic programming. Then we discuss 
the classical representation of CSP as a logic and con- 
straint logic program. The next section gives a different 
representation of a CSP and shows how it can be trans- 
lated to an equivalent logic program which then can be 
solved by abductive reasoning or by computing stable 
models of the program. Finally, we give the results of 
solving some well-known problems with several systems 
and the last section gives a summary of our research and 
some directions for future research. 

Preliminaries 

A constraint satisfaction problem is usually defined as 
a finite set of constraint variables X = {Xi, . . . , A„}, 
a finite domain Di of possible values for each vari- 



^We do not use the actual syntax of the systems but a 
syntax which is more uniform across the different systems. 



able Xi, and a finite set of constraint relations TZ = 
{rsi,--- ,rst}, where Si C {1,... ,n} for i = l,..,t 
are the indices of the variables which are constrained 
by rSi- A solution is an instantiation of the variables 
X which satisfies all the constraints in TZ. 

A logic program is defined for a fixed language C 
of function and predicate symbols and an infinite set 
of variables X. We assume that our language always 
contains the equality predicate = /2. Terms and atoms 
are defined as usual. A literal is either an atom a or a 
negation of an atom not a. A clause has the form 

a •*— 61, . . . 

where a is an atom and is called head and bi arc lit- 
erals and 61,... ,bn is called body. A clause is called 
definite if all bi are positive. A (definite) logic program 
P is a set of (definite) clauses. A Herbrand interpre- 
tation and Herbrand model are again defined as usual 
(cf. (Lloyd 1987)). A definite logic program always has 
a unique least Herbrand model which will be denoted 
with lm{P). 

Knowledge Representation of CSP 
Logic Programming 

A CSP can be easily represented as a definite logic pro- 
gram. The domains of the variables are defined by 
unary predicates. A constraint relation r/n is defined 

by a predicate p/n such that lm{p/n) = r/n. Then a 
CSP can be defined by a clause of the form: 

csp{Xi, . . . ,Xn) ^ 

di(Xi),... ,d„(X„), 
Pi{Xs^),... ,pt{Xs,). 

A computed answer substitution 6 for the goal <— 
csp{Xi, . . . , Xn) using SLD resolution will be a solu- 
tion to the CSP since P \= csp{Xi,... ,Xn)9 and in 
particular lm{P) \= csp{Xi,... ,Xn)6. It is also easy 
to see that the least model of the csp/n prc^dicate will 
contain all solutions to the problem. So. to compute 
all solutions to a given problem one can use a bottom- 
up fix-point operator to compute the least model of a 
program. 

To give an example let us consider the problem for 
finding the positions of n queens on a n x n chess board 
such that no two queens attack each other. One way 
of formalizing this problem is by using boolean vari- 
ables which indicate whether a queen is placed on a par- 
ticiilar square. However, by taking into account the fact 
that no two queens can be placed on the same column, 
we can use only n variables which give the positions of 
the queens for each column. Then for all pairs of queens 
we should state the constraint that they can not attack 
each other by being on a same row or diagonal. This 
constraint can be defined with the following predicate 
which is parameterized by the distance D between the 
columns of the two queens: 

sa/e(Xi, X2, D) ^ Xi 7^ X2, abs{Xi - X2) ^ D. 

Then the 4-Queens problem can be defined as: 



csp{Xi,X2,X^,Xi) 

d{Xi),d{X2),d{X3),d{X4), 

safe{Xi, X2, 1), saf e{X2,X3, 1), 
sa/e(Xi,X3,2),sa/e(Xi,X4,3), 
safe{X2,X4,2),safe{X3,X4,l). 

Here we exploit another particularity of the n-Queens 
problem, that the safe constraint is symmetric for the 
first two arguments. Hence it is enough to only check 
queens such that the column of the first one is less than 
the column of the second one. Executing the query 
<— csp{Xi,X2, X3, X4) under the top-down left- to-right 
strategy of Prolog will result in solving the problem by a 
generate and test approach. By interleaving the calls to 
the domain predicates d/1 and the constraint predicates 
will result in a standard backtracking: 

CSp{XuX2,X3,Xi) ^ 

d{Xi),d{X2), safe{Xi,X2A). 
d{X3), safe{X2, X3, 1), safe{X,, X3, 2), 
d{X4), safe{Xi, X4,3), safe{X2, X4, 2), 
safe{Xs,X4.l). 

A lot of research has been done to improve the ex- 
ecution strategy of standard Prolog. For example, a 
technique known as co-routing uses a literal selection 
strategy which selects a constraint as soon as all its ar- 
guments become ground. This allows a generate and 
test program to be executed with standard backtrack- 
ing. Another technique known as intelligent backtrack- 
ing (Bruynooghe 1991) docs a failure analysis and back- 
tracks only to variable assignments which are actually 
responsible for the failure. 

csp{N, L) <- 

make_vars{L, N), 
constrain jill[L) . 

mafce_uars([], 0). 
make_vars([H\T],N) 4- 

X > 0, 

d{H), 

Nl isN -1, 
make-vars{T, Nl). 

constrainjill{ [] ) . 
constrain-all{[X\X s\) <— 

constrain-all{Xs), 

constrain Jbetween{X , Xs, 1). 

constrain J}etween{X , [], X). 
constrain Jbetween{X , [Y\Ys],N) <— 
safeiX,Y,N), 

Nl isN + 1, 

constrainJ}etween{X, Ys, Ni). 

Figure 1: N-Queens as a Definite Logic Program 

An important problem with this representation of 

a CSP is that it has to be defined specifically for a 
given number of queens. If we want to parameterize the 



problem with respect to the number of queens then we 
should use some data structure (usually a list) to store 
the constraint variables. Figure 1 shows a typical spec- 
ification of the problem. The make-vars/2 predicate 
constructs a list of n variables and uses backtracking to 
enumerate all possible solutions. 

However this way of representing the problem is even 
less declarative than before. First of all, the column 
number of each queen is implicitly defined by the po- 
sition of the corresponding variable in the list. Then, 
in order to add constraints only between queens such 
that the column of the first one is less than the column 
of the second one we have to use two nested recursive 
predicates. 

graph-Coloring{Vert, Edges, Vars) <— 
make-vars{Vert, Vars), 
add-constr {Edges, Vars). 

make-vars{W, []). 

make-vars{[V\Vert], [assoc{V, X)\Vars]) <— 

color (X), 

makejjars{Vert, Vars). 
add-ccmstr{^, Vars). 

add-Constr{[edge{y\,V2)\Edges\,Vars) <— 

member [as soc{Vi, Xi), Vars), 
member{assoc{V2, X2), Vars), 
Xi X2, 

add.constr {Edges, Vars). 
Figure 2: Graph Coloring as a Definite Logic Program 

Another typical CSP problem is that of graph col- 
oring. The goal is to color all the vertices of a graph 
in such a way that no two adjacent vertices have the 
same color. The graph can be represented by a list 
of vertices and a list of terms edge{vi,V2) describing 
the edges. The problem can be expressed as a CSP 
by associating a different variable for each vertex and 
restricting its domain to all possible colors. Then for 
each edge in the graph we put a disequality constraint 
between the variables corresponding to the vertices of 
the edge. Figure 2 shows a sample formulation of this 
problem as a logic program. Here the correspondence 
between a vertex and the color assigned to it is made 
explicit by means of an association list. 

Constraint Logic Programming 

Constraint logic programming (CLP) (Jaffar & Maher 
1994) is an extension of logic programming where some 
of the predicate and function symbols have a fixed inter- 
pretation. This interpretation is dependent on a par- 
ticular constraint domain X (e.g. finite trees or real 
numbers) and this allows for a much more natural rep- 
resentation of problems from this particular domain. 
Besides its better declarative semantics for expressing 
problems, an implementation of a CLP(X) system also 



includes an efficient domain specific solver solvex for 
checking satisfiability of a set of constraints. 

A proof procedure for CLP is defined as an extension 
of standard resolution. A state is defined as a pair 
a, A II C) of a goal and a set of constraints. At each 
step of the computation, some literal a is selected from 
the current goal according to some selection function. 
If a is a constraint predicate then the next state is (<— 
^ II C A a) if solve{C A a) ^ false or (□ || false) 
otherwise. If a is a normal atom then the next state is 
s — t, B, A II C) for some clause b ^ B where a is 
of the form p{s) and b is of the form p{t). 

A well suited constraint domain for representing 
CSPs is that of finite domain integer arithmetic 
CLP{FD) (Codognet & Diaz 1996; Van Hentenryck, 
Saraswat, & Deville 1998). It includes standard arith- 
metic relations like =, 7^, < and fimctions +, — . * with 
their usual interpretation. The implementation of 
CLP{FD) is based on consistency (Tsang 1993) algo- 
rithms. 

A CSP represented as a logic program can be trans- 
lated to a CLP(FD) program in a straightforward way. 
First, the domains of the variables are declared by a 
special CLP predicate (for example X in l..n) and 
the constraints are defined using the CLP predicates 
and functions. Then the execution of the goal <— 
csp{Xi, . . . ,Xn) by a constraint proof procedure will 
result in a set of constraints which are then solved effi- 
ciently by consistency techniques. 

An Alternative Representation 

From the examples and the discussion in the previous 
section, the following general methodology can be given 
for representing a parameterized CSP as a (constraint) 
logic program. First, we create some data structure 
where each value of the parameter of the problem (e.g. 
the number of queens, or the number of vertices) is as- 
sociated with a different variable. Then we define some 
recursive predicates which iterate over this data struc- 
ture and define constraints between the variables whose 
corresponding parameters satisfy certain conditions. 

It seems natural to represent the mapping between 
parameter and variable with a function. First let us 
introduce a new domain Dp which includes all values of 
the parameter. For example, for the n-Queens problem 
Dp is the set of columns and for the graph coloring 
problem Dp will contain all vertices. Then a solution 
to a CSP can be represented by a function f : Dp ^ D, 
where D is the domain of the constraint variables. 

The use of functions for representing CSPs can be 
realized in a first-order logic setting where a solution 
of a problem can be given by an interpretation of a 
function symbol in a model of the theory describing 
the problem. Also the constraints are expressed more 
naturally by using function symbols. A definition of 
the n-Queens problem in first-order logic could be the 
following: 



VCi,C2. Ci < C2 ^ 

safe{pos{Ci),pos{C2),C2 - C\). 

VXi,X2,D. safe{Xi.X2,D) ^ 
Xi^X2,abs(Xi-X2)^D. 

The function symbol pos/1 represents the mapping 
from columns to rows. A modeP of this theory based 
on a domain with n elements will consist of an inter- 
pretation of the function pos/1 and will be a solution 
of the problem with n queens. 

If we want to work with several domains with differ- 
ent sizes then wc may use a many sorted first-order logic 
where the arguments of the function and predicate sym- 
bols are assigned (possibly different) sorts. Then in an 
interpretation of our theory for each sort we can assign 
a domain with the appropriate size. The full specifica- 
tion of the graph coloring problem in many sorted logic 
is given below. 

F = {col : Sv Sc} 

yVuV2. col{Vi) col{y2) ^ edge{VuV2). 
edge{l,2). 

Here we use two sorts - for vertices, and Sc for 
colors and a function col : ^ Sc which maps vertices 
to colors. The edges of the graph are described as a set 
of facts edge{vi,V2)- An interpretation of this theory 
should associate the set with all vertices of the graph 
with the sort Sy and a set with possible colors with the 
sort Sc- 

Several systems for generation of finite models of 
many sorted first-order logic theories are available (e.g. 

FINDER (Slancy 1995) and SEM (Zhang & Zhang 
1995)) and thus can be used for solving CSPs. 

Back to Logic Programming 

Using functions for representing CSPs in a logic pro- 
gramming setting is not possible. The reason is that 
the domain of the computation is the Herbrand do- 
main (or comes from the constraint domain) and that 
the interpretation of the function symbols is fixed. To 
overcome this restriction, one can introduce miary pred- 
icates dp/1 and d/1 defining the domain and the range 
of the functions and encode a Tinary function f{X) = Y 
as a binary predicate p j {X, Y) with domain defined by 
dp /I and range by d/1. The following axioms establish 
that the interpretation of Pf /I corresponds to a func- 
tion: 

VX. 3F. d{Y)hpf{X,Y) ^ dp{X). (1) 
\fX,Y. dp{X)^Pf{X,Y). (2) 
yX,Y,Z. Y = Z ^Pf{X,Y),pf{X,Z). (3) 

^We assume here that the arithmetic relations and func- 
tions have a fixed interpretation, lilce in CLP. 



The formula (1) states that the predicate pf{X,Y) is 
defined for all X in the domain dp{X) and that the value 
of Y is in the range d{Y). The formula (2) enforces 
that the predicate p/ is false for all values of X not in 
the domain dp(X) and (3) ensures that the predicate 
Pf{X, Y) defines a function from X to Y. By introduc- 
ing an auxiliary predicate hasjp, (1) can be rewritten 
as: 

VX. has4){X)'^3Y. d{Y) Apf{X,Y). (la) 
VX. has.p{X) ^ dp{X). (lb) 

Formula la is the completion of the predicate hasjp. 
Writing it as a logic program clause, the only-if part can 
be dropped. However, the other three axioms (lb), (2), 
and (3) do not define any predicates but are just formu- 
lae which should be true in models of the program. In 
logic programming, such clauses are known as integrity 
constraints and we will denote them with the symbol 
<=^. Figure 3 shows the resulting logic program. 

has.p{X) ^ d{Y),pf{X, Y). from (la) 

has4){X) <^ dp{X). from (lb) 

dp{X) ^Pf{X,Y). from (2) 

Y = Z ^Pf{X,Y),pf{X,Z). from (3) 

Figure 3: A logic program limiting a predicate to be a 
function. 

A solution to a CSP problem is then given by an 
interpretation of the predicate Pf/2 in a model of the 
logic program obtained by replacing the functions //I 
with the predicates Pf/2 and adding the theory defin- 
ing pf/2. Different models give rise to different so- 
lutions. Applying this transformation on the n-Qucens 
problem, one obtains the specification given in figure 4 
where pos{C, R) is the predicate which gives the row R 
of a queen at a column C. 

Abduction 

Abductive logic programming (Kakas, Kowalski, & 
Toni 1992) is a form of reasoning in which an answer to 

a query is a set of facts. More formally, an abductive 
framework is defined as a triple (P, A, I) where P is a 
logic program, A is a set of predicates called abducibles 
and J is a set of integrity constraints. A solution is a 
set A C such that M \= I for some canonical model 
M of P U A. If one is interested in a query, then it can 
be put as part of the integrity constraints. Different 
choices for the type of canonical model have been con- 
sidered in the literature. In (Kakas & Mancarella 1990), 

'^The standard form to write integrity constraints is as 
clauses with an empty head. Here we use another notation 
for better clarity. 

''Here and in the rest of the paper we will use the same 
symbol A to indicate both the set of abducible predicates 
and the set of all their ground instances. 



safe{Xi, X2, D)^Xi^ X2, abs{Xi - X2) ^ D. 
has.pos{X) <— drowO^)jPos{X,Y). 

has-pos(X) <= dcoi{X). 
dcoi{X)<=pos{X,Y). 
Y = Z^ pos{X, Y),pos{X, Z). 

safe{Ri,R2,C2-Ci) ^ 

pos{Ci,Ri),pos{C2,R2),Ci < C2. 

Figure 4: N-Queens as an (Abductive) Logic Program 



M must be a stable model of P U A which is also called 
generalized stable model M{A) of P. In (Denecker & 
Dc Schrcyc 1998) three valued models of the comple- 
tion of the program are considered and the abductive 
predicates must have a two- valued interpretation. 

A problem defined as in figure 4 can be given directly 
to an abductive procedure by declaring the predicate 
pos{X, Y) as abducible. 

Proof procedures for abduction are defined in a simi- 
lar way (Eshghi & Kowalski 1989) to, or as an extension 
(Denecker & De Schreye 1998) of SLDNF resolution. In 
each state of the derivation they also maintain a set A 
of already abduccd atoms. When an abductive atom is 
selected in the current goal, it is checked if it can be 
resolved with an atom from A. If this is not the case 
then the selected atom is added to the set of abducibles 
and a consistency derivation is started which checks the 
integrity constraints to see if this assumption does not 
lead to a contradiction. If the abducible atoms contain 
variables then during an abductive step these variables 
are replaced with skolem constants and the miification 
algorithm is extended to deal with them (Denecker & 
De Schreye 1998). We have already discussed the prob- 
lems of the SLD resolution for efficiently solving CSPs 
thus we can expect that the performance of such more 
complex proof procedures will be even worse. 

Recently, abduction has been extended to constraint 
logic programs. One of the main ideas is that the skolem 
constants which arc added as arguments of the abduced 
predicates are in fact existentially quantified constraint 
variables and their values can be computed by a con- 
straint solver. This allows us to have a more declarative 
representation of CSPs and still use efficient techniques 
for computing their solution. The first such integra- 
tion is the ACLP system of Kakas (Kakas & Michael 
1995) which is based on the proof procedure of (Eshghi 
& Kowalski 1989). Originally, ACLP was defined only 
for definite programs and integrity constraints and in 
(Kakas, Michael, & Mourlas 2000) it was extended to 
deal with negation as failure through abduction in a 
similar way as in (Eshghi & Kowalski 1989). A more 
recent integration of abduction and CLP is the SLD- 
NFAC system (Denecker & Van Nuffelen 1999) which is 



based on the abductive procedure SLDNFA (Denecker 
& De Schreye 1998). 

Logic Programming with Stable Model 
Semantics 

Another way of using "open" predicates in a logic pro- 
gram is to use the stable model semantics (Gelfond & 
Lifschitz 1988). A predicate p can be defined as having 
an open interpretation by the following two rules: 

p ■*— not p. 
p ■*— not p. 

where p is a new predicate symbol. This program has 
two stable models - {p} and {p}. If A is a set of pred- 
icates then we define with T{A) the following set of 
clauses 

T{A) = { p{x) ^ not p{x). \ pe A}U 
{ p{x) ^ not p{x). \ p £ A} 

Then to compute a stable model of a CSP as repre- 
sented in section 4, one only needs to add to the pro- 
gram the clauses T{A), where A is the set of predicates 
which are the result of the transformation of the func- 
tion symbols"''. 

In fact, there is a very strong relation between the se- 
mantics of an abductive framework and the stable mod- 
els of an equivalent logic program. It has been shown 
in (Satoh & Iwayama 1991) that M(A) is a general- 
ized stable model of an abductive framework (P, A^ I) 
iff there exists a stable model M' of (P U T{A),I) such 
that M' = M(A)UV where V = { p{x) \ p{x) e A\A}. 

The rules declaring the open predicates can be com- 
bined with some of the integrity constraints to obtain 
a more compact representation of the problem. For ex- 
ample, the integrity constraint stating that the open 
predicate pf should be false for all values not in the 
domain dp 

dp{X)^Pf{X,Y) 

can be omitted by using the rule 

Pf{X,Y)^dp{X),notpf{X,Y). 

The reason is that in the stable model semantics an 
atom can be true only if it is a head of some clause. 
The full specification of the n-Queens problem is given 
in figure 5. 

An efficient implementation of the stable model se- 
mantics is the smodels system (Niemela & Simons 
1996). It works with propositional rules and a special 
pre-processing program is used for grounding function- 
free range-restricted logic programs. The implemen- 
tation of the system combines bottom-up inference 
with backtracking search and employs powerful pruning 
methods. 

''In fact, exactly the same methodology for representing 
CSPs is proposed in (Niemela 1999) 



pos{X, Y) ^ dcoi{X), not pos{X,Y). 
pos{X, Y) ^ not pos{X, Y). 

safe{Xi,X2,D) Xi ^ X2,abs{Xi - X2) ^ D. 
has-pos{X) <— drow{Y),pos{X,Y). 

has_pos{X) <^= dcoi{X). 

Y = Z^ pos{X, Y),pos{X, Z). 

safe{Ri,R2,C2-Ci)^ 

pos{Ci,Ri),pos{C2,R2),Ci < C2. 

Figure 5: N-Queens with Stable Models Semantics 
Other Formalisms 

In (Cadoli & Palopoli 1998) is presented a language 
called DATALOG*^^^*^ which is an extension of data- 
log where only some of the predicates are minimized 
and the interpretation of the others is loft open. The 
semantics of the language originates from the nonmono- 
tonic formalism of circumscription and is defined as the 
minimal Herbrand model of the program w.r.t. a fixed 
interpretation of the open predicates. It is proven in 
(Cadoli & Palopoli 1998) that the data complexity of 
deciding whether a query is not entailed by the pro- 
gram is NP-complete which means that one can express 
any CSP in this formalism. However, as the language 
does not contain negation, one can not use directly the 
methodology discussed above for representing CSPs. 

In (Cadoli et al. 1999) is defined the language NP- 
SPEC which is an extension of DATALOG*-^^^*^ and allows 
a more natural representation of problems. The main 
difference is that it supports special meta-declarations, 
called tailoring predicates, restricting the domain and 
the interpretation of the open predicates. The most 
simple one is of the form Subset{domain/n,pred) 
which defines pred/n to be an open predicate and 
its interpretation should be a subset of the interpre- 
tation of the predicate domain/ n. Another declara- 
tion is Partition{domain/n,pred, m) which states that 
the predicate pred/n should be partitioned in m sets 
which is, in fact, equivalent to a function with domain 
domain/n and a range 0..m — 1. The two other tailor- 
ing predicates arc Permutation{domain/n,pred) and 
IntFunc{domain/ n, pred, min..max) which express re- 
spectively that pred/n is a bijection from domain/n 
to 0.. domain/n — 1 and pred/n is a function with a 
range min..max. Another extension of NP-SPEC is the 
support of a predefined arithmetic functions and pred- 
icates. 

The formulation of the graph coloring problem in the 
NP-SPEC language is given below. The input of the 
program consists of facts node /I and edge/ 2 describing 
the graph. 

Partition {node/1, color, 4). 
<^ edge{Vi,V2), color{Vi,Ci), color{V2,C2). 



Experiments 

The Systems 

The finite domain CLP package is the one provided with 
SICStus version 3.7. Given the reputation of SICS- 
tus and of finite domain CLP. one can assume it offers 
state of the art technology for CSP solving and it is a 
good yardstick to judge the performance of other sys- 
tems. The abductive system ACLP (Kakas, Michael, 
& Mourlas 2000) is a meta interpreter written in Pro- 
log, runs on Eclipse version 4.2 and makes use of its fi- 
nite domain package. The abductive system SLDNFAC 
(Deneckcr & Van Nuffclen 1999) is also a meta inter- 
preter written in Prolog but runs on SICStus version 
3.7 and makes use of the SICStus finite domain pack- 
age. The model generator SEM (Zhang & Zhang 1995) 
version 1.7 is a fine tuned package written in C. smodels 
(Niemela & Simons 1996) version 2.25, the system for 
computing stable models is implemented in C-|— I- and 
the associated program used grounding is Iparse version 
0.99.48. All experiments have been done on the same 
hardware, namely Pentium II. 

All systems based on a finite domain constraint solver 
used a labeling strategy which first selects variables 
with the smallest domain and then the ones which par- 
ticipate in the highest number of constraints®. 

N- Queens 

Figure 6 gives the running times for the different sys- 
tems and figure 7 gives the number of backtracks. The 
two abductive systems (ACLP and SLDNFAC) do not 
introduce any extra choice points compared to CLP and 
hence are not plotted in figure 7. Not surprisingly, CLP 
gives the best results. SLDNFAC is second best and, 
although meta-interpretation overhead increases with 
problem size, deteriorates very slowly. SEM is third 
but runs out of memory for large problems (it needs 
about 120MB for 27 queens). This is probably caused 
by a not very good techniques for grounding the prob- 
lem and exploring the search space. The times given 
for SEM do not include time spend by the operating 
system in managing the memory which becomes con- 
siderable for the larger instances of the problem. ACLP 
performs substantially worse than SLDNFAC and de- 
grades more quickly for the larger problems. It can 
likely be attributed to the more experimental nature 
of the implementation, smodels performs very poorly 
on this problem, in particular when compared with its 
performance on the graph coloring problem. As can be 
seen from figure 7 the main reason seems to be the large 
number of backtracks it does. 

The CLP consistency techniques seem to be much 
less sensitive to the domains size, and this carries over 
to the abductive systems which reduce the problem to 
a CLP problem and then use the CLP solver to search 
for the solution. 



^This strategy is sometimes abbreviated to ifc. 
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Figure 7: N-Queens: Number of backtracks. 




Figure 8: Graph Coloring 



Graph Coloring 

We used a graph generator program which is 
available from address http://web.cs.ualberta.ca/ 

~joe/Coloring/Generators/generate .html. Wc did 
our experiments with planar undirected graphs which 
are known to be 4-colorable. The graphs were generated 
using a 20% probability of assigning arcs. This results 
in dense graphs with a stable behavior. For this prob- 
lem, the domain of the solution variables (the number 
of colors) remained the same and we have modified only 
the parameter of the problem (the number of vertices) 
and consequently the number of constraints (arcs). Fig- 
ure 8 gives the results of solving the problem with the 
different systems. Both axes arc plotted in a logarith- 
mic scale. On the x-axis we have put the number of 
arcs (constraints) instead of the number of vertices. 

Not surprisingly, CLP is the fastest system, smod- 
els is second best on this problem. We assume it is in 
part because of the very concise formulation. Using the 
so called technique of rules with exceptions (Niemela 
1999), the two rules needed to describe the space of 
candidate solutions also encode the constraint that the 
color is a function of the vertex. Hence there is only 
one other rule, namely the constraint that two adja- 
cent vertices must have a different color. The differ- 
ence with CLP is almost two orders of magnitude for 
the largest problems. SLDNFAC is slightly worse than 
smodels. Although meta-interpretation overhead tends 
to increase with problems size, the difference with smod- 
els grows very slowly. The model generator SEM de- 
teriorates much faster and runs out of memory for the 



larger problems. The fact that it grounds the whole the- 
ory is a likely explanation. The difference with smod- 
els supports the claim that smodels has better tech- 
niques for grounding. ACLP performs substantially 
worse than SLDNFA and also deteriorates faster. 

Conclusion 

The examples which we have considered in this paper 
are by no ways representative. However we think that 
they still show some interesting features and limitations 
of the considered systems. 

Consistency algorithms are a very efficient way for 
solving CSP. Constraint logic programming allows this 
techniques to be integrated in a natural and clear way 
to logic programs. However, as argued parameterized 
CSPs can not be represented in a declarative way as 
CLP(FD) programs. Using abduction or stable mod- 
els as the basis for logic programming allows the prob- 
lems to be represented in a more declarative way and 
the recent integration of abduction with CLP allows 
the same consistency techniques to be used for solving 
the problems. At the moment, such systems arc imple- 
mented as meta-interpreters on top of Prolog and they 
essentially reduce a problem to the same set of con- 
straints (in many cases without backtracking) which 
would be produced by the corresponding constraint 
logic program. Our experiments suggest that the over- 
head of an abductive system is small and acceptable. 
Moreover they can also solve other classes of problems 
which require non-monotonic reasoning like planning 
problems. We also showed that there is a very close 



relation between the semantics and the problem rep- 
resentation of abduction and logic programming with 
stable model semantics. The only difference is in the 
reasoning techniques - abduction is usually done by 
a top-down proof procedure, while a stable model is 
usually computed by a bottom-up procedure. How- 
ever, the techniques used to compute a stable model 
of a program do not seem to be so well suited for 
solving CSPs. One reason could be that they work 
on a ground propositional programs which tend to be 
large and grow fast as the parameter of the problem in- 
creases. This suggests that an interesting area for fur- 
ther research would be a framework for computing sta- 
ble models of constraint logic programs with the help of 
constraint solving techniques. Some work has already 
been done in this direction (Dix & Stolzenburg 1998; 
Eiter, Lu, & Subrahmanian 1997). 

We argued earlier in the paper that the most nat- 
ural way for representing CSPs is with functions with 
open interpretation. Hence it would be interesting to 
consider extensions of the CLP scheme which directly 
support such form of reasoning. Some work in this 
area is done in (Bruynooghe, Pelov, & Denecker 1999; 
Rickey 1993). 
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